Four-point vertex in the Hubbard model and partial bosonization 
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Magnetic and superconducting instabilities in the two-dimensional -Hubbard model are discussed within a 
functional renormalization group approach. The fermionic four-point vertex is efficiently parametrized by means 
of partial bosonization. The exchange of composite bosons in the magnetic, charge density and superconducting 
channels accounts for the increase in the eff'ective couplings with increasing length scale. We compute the 
pseudocritical temperature for the onset of local order in various channels. 
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I. INTRODUCTION 

The functional renormalization group approach to corre- 
lated fermion systems has been of great help to detect different 
types of instabilities and collective order within many differ- 
ent models. This holds, in particular, for the two-dimensional 
Hubbard model which is hoped to improve our understand- 
ing of superconductivity in the high-Jc-cuprates.! 1-8] Most 
studies presented so far rely on the flow of the momentum- 
dependent four-fermion vertex. (For analogous work on four- 
quark vertices see f9', To^.) They are performed in the so- 
called A/^-patch scheme where the Fermi surface is discretized 
into N patches, and the angular dependence of the fermionic 
four-point function is evaluated only for one momentum on 
each directional patch. 

The approach presented here is based on the introduction 
of fermionic bilinears corresponding to different types of pos- 
sible orders (partial bosonization) that was developed and 
used before in Illl4l6i1 . It is also inspired by the efficient 
parametrization method for the fermionic four-point vertex 
proposed in iflTh . The link between the two approaches is 
given by the fact that different channels of the fermionic four- 
point function, defined by their (almost) singular momen- 
tum structure, correspond to different types of possible orders 
which are described by different composite bosonic fields. 

The advantages of our method are, first, that it allows one 
to treat the complex momentum dependence of the fermionic 
four-point function in an efficient, simplified way, involving 
only a small number of coupled flow equations and second, 
that it permits to follow the renormalization group flow into 
phases of broken symmetry. A comparative disadvantage may 
be a better resolution of contributions from many channels in 
the A/^-patch approach. (In principle, both approaches can be 
combined.) In this paper we focus on the first of these two as- 
pects. Spontaneous symmetry breaking was already addressed 
for antiferromagnetism (AF) llj, [iSO in the Hubbard model 
close to half filling and for J-wave superconductivity in an 
effective Hubbard-type model with a dominating coupling in 
the J- wave channel. 1 14] 

The main idea behind partial bosonization, namely, to rep- 
resent the fermionic four-point vertex by a certain number of 
exchange bosons, is graphically shown in Fig. [H The de- 
pendence of the four-point vertex on three external frequen- 
cies and momenta (or simply "momenta", as we are going to 
write for short) is parametrized in terms of bosonic propaga- 



FIG. 1 : Schematic picture of bosonization of the four-fermion ver- 
tex. Solid lines correspond to fermions, the dashed line to a complex 
(Cooper pair) boson, the wiggly line to a real boson representing a 
particle-hole state in the spin or charge density wave channel. 



tors together with Yukawa couplings which describe the in- 
teraction between one boson and two fermions. The different 
bosonic channels are distinguished according to the structure 
of their momentum dependence which may possibly become 
singular due to a zero of the inverse boson propagator. The 
momentum-independent part of the four-fermion vertex may 
either be distributed onto the different bosonic channels or be 
kept fixed as a purely fermionic coupling. In order to avoid the 
arbitrariness encountered when one chooses the first of these 
two options, we adopt the second for our computations. This 
may be regarded as a prototype for a combination of partial 
bosonization with the A/^-patch method in the sense of keep- 
ing only one patch and setting the four-fermion coupling to a 
constant value. 

Although for a numerically exact treatment of the four- 
fermion vertex an infinity of bosonic fields would have to 
be considered in principle, a small number of well-chosen 
fields may suffice for a reasonable quantitative precision. The 
choice of fields that has to be included in order to capture 
the relevant physics depends on the model under investiga- 
tion. In the case of the two-dimensional Hubbard model at 
small next-to-nearest-neighbor hopping \ f\, a magnetic boson 
m and a J- wave Cooper pair boson d are needed because they 
correspond to the instabilities that occur. In order to avoid a 
too poor momentum resolution of the four-fermion vertex, we 
also include an ^-wave Cooper pair boson field s and a charge 
density boson field p. Other types of bosons are needed in 
other contexts — for instance, a J- wave charge density boson 
for the study of Pomeranchuk instabilities or a /7-wave boson 
for triplett superconductivity at larger values of |^'| away from 
the van Hove filling |7]. With the restriction to the bosons 
m, J, s and p, supplemented by a pointlike four-fermion ver- 
tex, we will show that interesting results and a semiquantita- 
tive understanding can be found in a rather simple truncation. 

The two-dimensional Hubbard model |[T8| - [2Q|] on a square 
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lattice has attracted a lot of attention in the past 25 years be- 
cause it is thought to cover important aspects of the physics of 
the high-Tc cuprates. In analogy to the phase diagram of the 
cuprates, it shows antiferromagnetic order at half-filling and is 
believed to exhibit J- wave superconducting order (dSC) away 
from half filling. 1I21I] Today there are many studies which pre- 
dict the J- wave instability to be the dominating one in a cer- 
tain range of parameters aside from half filling [22- 35], for a 
systematic overview see |36|] . The picture is also confirmed by 
some strikingly simple scaling approaches ll37l-[391 and finds 
further support within more elaborate renormalization group 
studies such as 



II. METHOD AND APPROXIMATION 

The starting point of our treatment is the exact flow equa- 
tion for the efl'ective average action or flowing action, 

dkTk = ^STr (rf + Rk)"' ^kRk = ^STr 4 (ln(rf + Rk)) . 

(1) 

The dependence on the renormalization group scale k is intro- 
duced by adding a regulator Rj^ to the full inverse propagator 
Tf\ In Eq. ([T]), STr denotes a supertrace, which sums over 
momenta, frequencies, and internal indices, while is the 
scale derivative acting only on the infrared (IR) regulator Rj^. 
The Hamiltonian of the system under considerations is taken 
into account by the initial condition Tk^A = S of the renormal- 
ization flow, where A denotes some very large UV scale and S 
is the microscopic action in a functional integral formulation 
of the Hubbard model. In the IR limit (k 0) the flowing ac- 
tion Tk equals the full efl'ective action T = Tk^o, which is the 
generating functional of one-particle-irreducible (IPI) vertex 
functions. 

We employ a compact notation with Q = (ojn = 27rnT, q) 
and Q = (oJn = (2n -h l)7iT, q) for bosonic and fermionic fields 
and 

Q n--oo ^ V / 

m - 20 = r-i^,,,.(27r)V2\q - qO . (2) 

The components of the momentum q are measured in units of 
the inverse lattice distance a"^ The discreteness of the lattice 
is reflected by the 27r-periodicity of the momenta q. 

Although Eq. ([T]) is an exact flow equation, it can only 
be solved approximately. In particular, a truncation has to 
be specified for the flowing action, indicating which of the 
(infinitely many) IPI vertex functions are actually taken into 
account. Our ansatz for the flowing action includes contribu- 
tions for the electrons, for the bosons in the magnetic, charge, 
and ^-wave and J- wave superconducting channels, and for in- 
teractions between fermions and bosons, 

^k\x\ = ^F,klx] + ^Fm,klx] + ^Fp,klx] + ^Fs,klx] + ^Fd,kM 

^r^,klx] + r^M + r^M + ^d,klx] • (3) 



The collective field ;^ = (m,p, s, s*, d, J*, i//, il/^) includes both 
fermion fields ij/^ and boson fields m,p, s, s* , d, J*. 

The purely fermionic part Tflx] (the dependence on the 
scale k is always implicit in what follows) of the flow- 
ing action consists of a two-fermion kinetic term F/rkin, 
a momentum-independent four-fermion term F^, and the 
momentum-dependent four-fermion terms F^, F^, F^ and F^, 

Tf M = r^n + F^ + F^ + F^ + F^ + F^ . (4) 
The fermionic kinetic term is given by 

F^km=X ^^(0/^^(0^(2), (5) 

Q 

with inverse fermion propagator 

PF(Q) = ioj^m, (6) 

where we take for the dispersion relation of the free electrons 

^(q) = -jd- 2^(cos qx + cos qy) - 4/ cos qx cos qy . (7) 

The momentum-independent part of the four-fermion cou- 
pling is identical to the Hubbard interaction U. In our trunca- 
tion, this coupling is not modified during the flow. The corre- 
sponding part of the efl'ective action therefore reads 

^F = \ Yj US(Ki-K2 + K3-K4) 

X [ilf\Ki)ilf{K2)] [ilf\K3)ilf{K,)] . (8) 

In this work, as in contributions to the fermionic 

self-energy are neglected. Instead, we focus on the mo- 
mentum dependence of the fermionic four-point function 
Af(Ki,K2, ^3, ^4), which, due to energy-momentum conser- 
vation K4 = Ki - K2 -\- K3, is 3. function of three indepen- 
dent momenta. We decompose this vertex into a sum of four 
functions /l^(2), /l^(2), A'p(Q) and Af(Q), each depending on 
only one particular combination of the Ki, which correspond 
to the four difl'erent bosons taken into account. This is inspired 
by the singular momentum structure of the leading contribu- 
tions to the four-fermion vertex. In our ansatz for the efl'ective 
average action these functions enter as 

F^ = Yj "^FiKi - K2)S(K, -K2 + K,- K,) 

X {ily\K,)<7,lf{K2)\ ■ [,lf\K^)a,lf{K^)\ , (9) 

^F = -\ Z " S {Ki -K2 + K3- K4) 

K\ ,K2,K2i,K4 

X [,{^HKiMK2)] [^^(^3^(^4)1 (10) 

for the real bosons, and, for the superconducting bosons, as 
= Yj "^f^^' +Ki)5{Ki-K2 + K3- K4) 

X [^\Ki)eriK3)] [^^(/s:2)e<A(^4)] , d D 
r^= Yj 4(Ki+K3)6(Kr-K2 + K^-K4) 

K\ ,K2,Kt,,K^ 

xM{Ki-K,)/2)M{K2-K4)/2) 
X [^^.+(^1)6^(^3)] [iJf''(K2)e,Jf(K4)] , (12) 
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where cr = (cr^ , cr^, cr^) is the vector of the PauH matrices, the 
matrix e is defined as 6 = icr^, and the function 

fd(Q) = fd(q) = ^ (cos^;c - cos qy) (13) 

is the J- wave form factor which is kept fixed during the flow. 

In a first step, contributions to the four-fermion vertex are 
distributed onto the coupHngs A^, /l^, /l^, A^, depending 
on their momentum dependence. Partial bosonization comes 
into play at this stage as the absorption of these contribu- 
tions by the corresponding Yukawa couplings and bosonic 
propagators. More concretely, this means that the couplings 
Ap, /l^, A^j., Af. are set to zero by introducing a scale- 
dependence of the bosonic fields, which in turn generates ad- 
ditional contributions to the various Yukawa couplings. The 
technique by means of which this is achieved is called flow- 
ing bosonization or rebosonization.|43, 44] We describe it in 
some detail in Appendix A. In consequence, the complicated 
spin and momentum dependence of the fermionic four-point 
function Af(Ki , K2, K^, K4), Sisit emerges during the flow, will 
be captured by the momentum dependence of the propagators 
of the bosons and the couplings between bosons and fermions. 

The interaction between electrons and composite bosons 
are taken into account in our ansatz for the flowing action by 
Yukawa-type vertices of the form 

TFm = -Yj hm(K)m(K)-[ilf\Q)cTilf(Q')]S(K-Q+Q'), 

K,Q,Q' 

^Fs = -Yj hsiK) [s\K)\ilf\Q)eilf{Q)^ (14) 
-s(K) [i/r\Q)er(&)]) S(K-Q- Q) , 

= w^)/a(G-eo/2)(j*(/r)[iA^(e)6iA(eo] 

-d{K) V^\f\Q)erm^ S{K-Q- Q') . 

Note the presence of the J-wave form factor in the second- 
to-last line. To determine the ^-dependence of the Yukawa 
couplings hm, hp, hs, hd is a central task within our approach. 

The purely bosonic parts of the efl'ective action are char- 
acterized by the bosonic propagators. For the magnetic bo- 
son, for instance, the inverse propagator is given by Pm(Q) = 
Pm(Q) where is its minimal value and Pm(Q) is the 

(strictly positive) so-called kinetic term. The contributions to 
the efl'ective average action where the bosonic propagators ap- 
pear are 

= i^m^(-e)(p^(e) + m2^m(e), (15) 

^ Q 

= ^J]^(-Q){Pp(Q)^ml)p(Q), (16) 
^ Q 

r. = ^^*(G)(/^.(G) + m?)^(e), (17) 

Q 

Td = Yj'^\Q){Pd(Q)^m])d(Q). (18) 
Q 



Our parametrization of the frequency and momentum depen- 
dence of the bosonic propagators and the Yukawa couplings is 
described in Appendix B. In contrast to the decomposition of 
the fermionic four-point vertex proposed in [17], our bosonic 
propagators exhibit an explicit frequency dependence. 

In the present paper, the purely bosonic parts of the flowing 
action are confined to the bosonic propagators. Higher order 
purely bosonic interactions are currently investigated and will 
be included in a forthcoming work. 

III. INITIAL CONDITIONS AND REGULATORS 

At the microscopic scale ^ = A the flowing action must be 
equivalent to the microscopic action of the Hubbard model, so 
the initial value of the four-fermion coupling must correspond 
to the Hubbard interaction U. The bosonic fields decouple 
completely at this scale, so the initial values of the Yukawa 
couplings are 

hmk = hp\A = hs\A = hd\A = . (19) 

The purely bosonic part of the eff'ective action on initial scale 
is set to 

r^|A = m^.m, TplA =p^p, (20) 
r,|A = s*s, TdlA = d^'d. 

In other words, we take frtf^ = f and then use units t = I and 
Pi^A = 0. The choice m^^ = amounts to an arbitrary choice 
for the normalization of the bosonic fields, which are intro- 
duced as redundant auxiliary fields at the scale A, where they 
do not couple to the electrons. Of course, this changes dur- 
ing the flow, where the bosons are transformed into dynami- 
cal composite degrees of freedom, with nonzero Yukawa cou- 
plings and a nontrivial momentum dependence of their prop- 
agators. 

In addition to the truncation of the efl'ective average ac- 
tion, regulator functions for both fermions and bosons have 
to be specified. We use "optimized cutoff's" |45, 46] for both 
fermions and bosons. The regulator function for fermions is 
given by 

Rk (Q) = sgn(^(q)) (k - l^(q)l) 0(k - |^(q)|) , (21) 
the regulator functions for the real bosons are given by 

K^'iQ) = A^ip . {k^ - F,/,-(q, q))m^ - Fc//(q, q)) (22) 

allowing for an incommensurability q with Fcn as defined in 
Appendix B. Regulator functions for the Cooper-pair bosons 
are of the same form, but no incommensurability needs to be 
accounted for in these cases. 



IV. FLOW EQUATIONS 

The ffow equations for the couplings follow from projec- 
tion of the ffow equation for the ffowing action onto the var- 
ious diff'erent monomials of ffelds. The right-hand sides of 



O*" 

FIG. 2: IPI diagrams contributing to the flow of bosonic propa- 
gators. Wiggly lines denote real bosons (particle-hole channels), 
dashed lines complex bosons (Cooper pair channels). 




FIG. 3: IPI diagrams which directly contribute to the flow of the 
Yukawa couplings. 



these flow equations are given by the one-particle-irreducible 
diagrams having an appropriate number of external lines, in- 
cluding a scale derivative dk acting only on the IR regulator 
Rk. Diagrams contributing to the flow of boson propagators 
are shown in Fig. [21 

Once some bosonic mass term changes sign from posi- 
tive to negative during the flow, this signals the divergence of 
the four-fermion vertex function in the corresponding chan- 
nel. A negative mass term indicates local order, since at a 
given coarse graining scale k the efl'ective average action eval- 
uated at constant field has a minimum for a nonzero value of 
the boson field. The largest temperature where at fixed values 
of U,t\/j one of the mass terms fhf changes sign during the 
flow is called the pseudocritical temperature Tpc. It can also be 
described as the largest temperature where short-range order 
sets it. If this order persists for k reaching a macroscopic scale, 
the model exhibits efl'ectively spontaneous symmetry break- 
ing, associated in our model to (either commensurate or in- 
commensurate) antiferromagnetism or J-wave superconduc- 
tivity. The largest temperature for which local order persists 
up to some k corresponding to the inverse size of a macro- 
scopic sample is the true critical temperature T^. In this paper 
we focus on the symmetric regime where we have a positive 
mass term and stop the flow once a mass term reaches zero. 
We plan to address the symmetry-broken regimes in a future 
work. 

The flow equations for the Yukawa couplings consist of 
a direct contribution and an "indirect" contribution resulting 
from flowing bosonization, see Appendix A. Diagrams con- 
tributing directly to the flow of the Yukawa couplings are 
shown in Fig. [3] those that contribute via flowing bosonization 
are displayed in Fig. (H Since we choose to distribute contri- 
butions from flowing bosonization only onto the Yukawa cou- 
plings and not onto the masses, it is crucial to include a mo- 
mentum dependence of the Yukawa coupling hm in the mag- 
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FIG. 4: IPI diagrams contributing to the flow of the Yukawa cou- 
plings via flowing bosonization. 

netic channel in order to account for the emergence of the d- 
wave superconducting instability. Otherwise the contribution 
of the particle-particle box diagram (the first in the lower line 
of Fig. 131) to the J- wave coupling would be underestimated. 

In order to demonstrate how the contributions to the four- 
fermion vertex are taken into account via flowing bosoniza- 
tion, we discuss the case of the purely fermionic loop dia- 
grams shown in the upper line of Fig. |4l As long as no scale 
dependence due to the regulator function has been introduced, 
they are given by 

Ki,K2,K3,K4 P 



Pf(P)Pf(P + K3) Pf(P)Pf(-P + /^i + K^)! 
S (Ki -K2^K3- K4) [i/i\K,)ilf(K2)] • [i/iHK,)i/,(K,)] . 

In order to obtain the resulting contribution to the fermionic 
four-point vertex function Ar^^^\ we have to take the fourth 
functional derivative of AF^ with respect to the fermionic 
fields. It is given by 

AY^/'\KuK2,K3,K,)= (24) 



4Sra(Ki)Si/f^(K2)Sry(K3)Si/fs(K4) 

- "T 2^\pf(p)Pf(-p + + K3) 

~ Pf{P)Pf{P + K2-Ki) ^ Pf(P)Pf(P^K2-K3V ' 

where S ay;/36 = \ [^a/3^y6 " ^ad^y/s) dcnotcs the singlct projec- 
tion. The two last lines of Eq. ((241) can be compared to the 
fourth derivative with respect to the fields of the right hand 
sides of Eqs. ^ - ([T2l) . This allows one to obtain the loop 
corrections to the four-fermion couplings A^, /l^, /l^, 
introduced there. The second last line of Eq. ((241) can be ab- 
sorbed by the ^-boson, the last line by the m- and p-bosons. 
No contribution to the J-boson arises at this stage. 
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FIG. 5: Schematic picture of the bosonization of loop contributions 
to the four-fermion vertex. The terms indicated by the three dots 
correspond to loop diagrams having internal bosonic lines. 



To determine how the last line of Eq. (l24l) should be dis- 
tributed onto the m- and p-bosonic channels, we use the iden- 
tity SadSyp = \ {SafiSys + ^ap^yd) ' tcrms havc now the 
same structures as those appearing in the fourth functional 
derivative of Eq. We obtain the following loop contri- 
butions to A"^, 4, 4: 



Y 1 

'2 ^ Pf(P)Pf(P^K2-Ki) 



2 ^ Pf(P)Pf(P^K2-Ki) 



,(25) 



2 ^ Pf(P)Pf(-P^Ki^K,) 



Ap, 



Ap is obtained from the one- 



The k dependence of A^, 

loop expressions (l25l) by adding the infrared cutoff to 
the inverse fermionic propagator and by applying the formal 
derivative dk = (dkR^)d/dR^ under the summation. For A^, 
for example, one obtains 



dkA'^(Q) = dkAA'^(Q), 



(26) 



where the formal derivative dk should be read as acting under 
the loop summation of terms contributing to AAp(Q). Note 
that (AAp)^(Q) is only part of the complete loop contribution 
AA^(Q), namely, the one which arises from the two diagrams 
shown in the first line of Fig. IH The complete AAp(Q) is 
obtained if the diagrams shown in Fig. IHare all taken together. 

In our partially bosonized approach, the fermion loop con- 
tributions to the momentum-dependent four-fermion vertex in 
Eq. (l25l) are fully accounted for by the exchange of the bosons 
m,p and s. This is shown schematically in Fig. \5\ In the lan- 
guage of boson exchange, the momentum dependence of the 
coupling in, for instance, the magnetic channel can be taken 
into account by the momentum dependence of the expression 
hl(Ki - K2)P-^{Ki - K2). In practice, we keep = during 
the flow and account for the loop-generated AA^ by a cor- 
responding change Ah^. We note that only the combination 
h^P'^ appears in the computations as long as the only mo- 
mentum dependence of the Yukawa couplings is that of the 
boson momentum. In fact, by a momentum-dependent rescal- 
ing of the fields for the m-boson it is in principle possible to 
arbitrarily attribute parts of the momentum dependence to 
or to Pm- Nevertheless, introducing the two factors and P'^ 
instead of only A^ is useful if one wants to approach sponta- 
neous symmetry breaking. It has the advantage that instead of 
having to deal with the divergent coupling A^, one only needs 



to account for the mass term changing its sign. The term con- 
taining in Eq. (O, which is quadratic in the boson field, 
becomes part of the eff'ective potential for the magnetic bo- 
son in a more extended truncation. Our description of the 
^-dependent flow of A^ by means of the ^-dependent quan- 
tities hm and Pm (and analogously for 4 is achieved 
formally by a /:-dependent nonlinear field redefinition, see Ap- 
pendix A, Eq. (IA5b . At momentum 2 = 0, for example, the 
contribution to the flow of the momentum-dependent Yukawa 
couplings due to the diagrams in the first line of Fig. |4l ac- 
cording to Eqs. (IA7b . (IA8b , is given by 



Pf(P)Pf(-P) 



(28) 



At this level, we have described the exact one-loop pertur- 
bative result for the momentum-dependent four-fermion ver- 
tex in terms of boson exchange. The concept of the flowing 
action, however, allows for a "renormalization group improve- 
ment" which is obtained by ^-dependent "running couplings" 
or vertices. In the purely fermionic flows [1- 8] the constant 
coupling U would be replaced by the full momentum- and k- 
dependent four-fermion vertex. In our partially bosonized ap- 
proach, where we keep only a constant four-fermion coupling 
U, this renormalization group improvement is generated by 
the diagrams involving internal bosonic lines, shown in Fig. [3] 
and the second and third lines of HI It is at this level where our 
truncation for the momentum dependence of the Yukawa cou- 
plings and inverse boson propagators as well as the restriction 
to a certain number of bosons starts to matter. 

The momentum dependence of the four-fermion vertex 
which is generated by boson exchange is much more com- 
plicated than the simple form (|25]) . We therefore have to de- 
cide how to distribute these contributions onto the difl'erent 
boson exchange channels. To this end, we adopt an approxi- 
mation where the momentum-dependence of the four-fermion 
couplings A^, A^p, A^, Af. can be identified with the depen- 
dence of the diagrams in Figs. [3] and H] on the so-called trans- 
fer momentum. This momentum is defined as the diff'erence 
between the momenta attached to the two fermionic propaga- 
tors in each diagram. Particle-hole diagrams are absorbed by 
the real bosons and particle-particle diagrams by the complex 
Cooper pair bosons. All diagrams are evaluated at external 
momenta L = (nT, tt, 0) and U = (nT, 0, tt) and transfer mo- 
menta = (0, 0, 0) and IT = (0, tt, tt). For small values of 
and l^'l, the (spatial parts of) momenta L and L' are close to the 
Fermi surface and the density of states is rather large there, so 
that this choice will capture the relevant physics for not too 
large |yu| and Where more than one combination of exter- 
nal momenta ±L and ±U is compatible with the condition that 
the transfer momentum is either or IT, we take the average 
over them. For the coupling in the J- wave channel, the evalu- 
ation of the contributing diagrams is discussed in more detail 
in the next section. 

While the contributions to the Yukawa couplings in Eqs. 
(l27l) and ([28]) are proportional to and therefore present al- 



6 



ready for large k, the diagrams shown in Fig. [3] and in the 
second and third Hnes of Fig. |4l start to have an influence on 
the flow of the Yukawa coupHngs only after nonzero Yukawa 
couplings have been generated due to Eqs. (1271) and (1281) in 
the first place. In perturbation theory, they would correspond 
to higher order eff'ects ~ and [/^. (Perturbatively, every 
Yukawa coupling counts as U.) The flow of the couplings in 
the magnetic and charge density channels starts to difl'er as 
soon as the diagrams shown in the first line of Fig. [3] become 
important. They contribute positively to the coupling in the 
magnetic channel but negatively to the couplings in the charge 
density and superconducting ^-wave channels. This explains 
why among the three Yukawa couplings hm,hp, hs the domi- 
nating one is hm, although in accordance with Eqs. (1771) and 
(1281) all three are generated with equal size at early stages of 
the flow. Due to the comparatively large Yukawa coupling 
hm the mass term is driven fastest toward zero by the dia- 
grams in Fig. [21 We can therefore understand why the charge 
density and ^-wave superconducting channels never become 
critical in the range of parameters investigated. 



V. COUPLING IN THE J- WAVE CHANNEL 

The generation of a coupling in the J- wave channel in the 
framework used here has already been discussed in an earlier 
work. 1 15] The J- wave Yukawa coupling arises during the 
renormalization flow due to the first diagram in the lower line 
of Fig. m which is the only particle-particle box graph. The 
coupling is extracted from contributions due to this graph by 
means of the prescription 

A4(l,r) = h^^Yff\L^L,-L,-L) (29) 

where the subscript s denotes the singlet and the superscript 
pp the particle-particle part of the four-point vertex. The mo- 
mentum vectors L and L are defined as in the previous sec- 
tion. For a motivation of this definition of the J- wave cou- 
pling see 1 1^ . The contribution from the particle-particle box 
diagram to the 5-- wave superconducting channel is obtained by 
adding, instead of subtracting, the two terms on the right-hand 
side of Eq. (|29l ). The s- and J- wave superconducting channels 
of the four-fermion coupling can be described as those parts of 
its singlet particle-particle contribution which are symmetric 
(5-- wave) and antisymmetric (J- wave) under a rotation by 90° 
of the outgoing electrons with respect to the incoming elec- 
trons. In our approximation, the first diagram in the second 
line of Fig. HI contributes only to the ^-wave channel. 

Once a coupling in the J-wave channel has been gener- 
ated through the particle-particle box diagram, it is further en- 
hanced due to the direct contribution shown as the first graph 
in the second line of Fig. [3l Since this graph, which is itself 
proportional to the Yukawa coupling hd, contributes positively 
to the flow of hd, it can lead to a growth of this coupling with- 
out bounds, i.e. lead to an an instability in the J- wave chan- 
nel. This instability will be the result of antiferromagnetic 



spin fluctuations (corresponding to the wiggly internal line of 
the diagram mentioned), so that our results finding a J- wave 
instability through this contribution support the idea, proposed 
and defended in |'22"-58^, that antiferromagnetic spin fluctua- 
tions are responsible for J- wave superconductivity in the two- 
dimensional Hubbard model (and maybe also in the cuprates 
insofar as the Hubbard model serves as a guide to the relevant 
cuprate physics). 

That the particle-particle graph in the second line of Fig. 
[3l is crucial for the emergence of a J- wave instability aris- 
ing from antiferromagnetic fluctuations is mirrored by the fact 
that this diagram has the same momentum structure as the 
BCS gap equation. In the presence of an interaction which 
in momentum space is maximal around the (;r, 7r)-points — a 
condition which is fulfilled when antiferromagnetic spin fluc- 
tuations dominate — the gap solving this equation exhibits d- 
wave symmetry. 



VI. NUMERICAL RESULTS 

We now come to the discussion of the numerical results 
we have obtained at small next-to-nearest-neighbor hopping 
\f\lt < 0.1. For values of \f\ and |yu| which are larger than 
those for which we show results in Fig. [8l there is, in addi- 
tion to the tendency to antiferromagnetism which is present 
already at large scales ^ » 0, a tendency toward ferromag- 
netism which becomes important at lower scales. Due to our 
simple parametrization of the inverse magnetic propagator Pm, 
see Eqs. (IB 21) - (IB4I) in Appendix B, and due to our choice of 
external momenta for evaluating the diagrams shown in Figs. 
[SlandlH we may overestimate magnetic fluctuations at larger 
values of |yu| and \ f\. Hence, we do not show any results for 
these larger values. 

In the range of parameters investigated, we find that either 
antiferromagnetism or J- wave superconductivity is the lead- 
ing instability. In agreement with previous findings, the cou- 
pling in the J- wave channel emerges due to antiferromagnetic 
fluctuations. In the parameter regime where this coupling is 
enhanced most strongly it competes with (and is driven by) 
the coupling in the incommensurate antiferromagnetic (iAF) 
channel which was studied in detail in Ref. 16, where the 
same framework was used as here. 

In the upper panels of Figs. [6land[7l the flow of the difl'er- 
ent channels of the fermionic four-point vertex is displayed at 
fixed fit = -0.01 and diff'erent values of U/t = 2.5,3, 3.5. 
The lower panels show the flow of the corresponding bosonic 
masses, which approximate unrenormalized inverse suscepti- 
bilities in channels which are close to critical. As expected, 
the antiferromagnetic coupling grows fastest and remains the 
dominant one for small to intermediate values of |yu|, for a rep- 
resentative case see Fig. [6l For ji/t < -0.28, however, the 
J- wave coupling diverges for higher temperatures than the an- 
tiferromagnetic coupling, for an example of this kind of sce- 
nario see Fig. [71 The couplings in the charge density wave 
and superconducting ^--wave channels are also considerably 
enhanced in both cases, and their influence is quantitatively 
important although they do not diverge. 
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FIG. 6: Upper panel: flow of the four-fermion vertex in the dif- 
ferent channels for U/t = 3, f/t = -0.01, fi/t = -0.12 and 
T/t = 0.188. The shorthands used in the legend are defined as 
^F,ac ^ hi{U)/P,n{Ul Af,, ^ hl{U)/P^{Ul Af,, ^ mO)/PAO) and 
Af4 = h^(^)/Pd(0). Lower panel: flow of the minima of the inverse 
bosonic propagators (bosonic mass terms). 
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FIG. 7: Same as Fig.[6]for ji/t = -0.32 and T/t = 0.109. In addition 
to the couplings defined in Fig[6]we plot Ap^ai = ^m(n- 0/^m(n- Q) 
(coupling in the incommensurate AF channel), where Q = (0, §, 0) 
with q the size of the incommensurability. For these parameters the 
coupling in the J- wave channel diverges first. In the magnetic chan- 
nel, incommensurate antiferromagnetic fluctuations (long-dashed 
lines) dominate over commensurate ones (short-dashed lines). 



In Fig. [8] the highest temperature at a given value of ji is 
plotted for which one of the boson masses drops to zero at 
some scale k, signaling the onset of local order on a typical 
length scale in the corresponding channel. These "pseu- 
docritical temperatures" Tpc are shown for f/t = -0.01 (up- 
per panel) and f/t = -0.1 (middle panel) and different val- 
ues of U. Pseudocritical temperatures for antiferromagnetism 
are higher by a factor ^ 3 than those presented in our last 
paper. 1 16] This is mainly due to the neglect of the fermionic 
wave function renormalization and of quartic bosonic cou- 
plings in the present paper. Both of these would suppress the 
growth of the four-fermion vertex and hence the emergence of 
local order. These contributions are omitted here for the sake 
of a simple and nevertheless systematic approach to the four- 
fermion vertex. They will be included in a forthcoming work 
where also the bosonic vertex functions that directly couple 
together the different types of bosons will be taken into con- 
sideration. We recall that often the true critical temperature 
Tc is found to be substantially smaller than the pseudocritical 
temperature r^c-lllllll 

For some pairs of parameters U and f there exists a range 
of values of jd where incommensurate antiferromagnetism has 
the largest pseudocritical temperature, for others there are no 
such values of yu, see Fig. [8l In the range of yu, where the tran- 
sition from (either commensurate or incommensurate) antifer- 



romagnetic to J- wave superconducting order occurs in Fig. [SI 
there is an extremely close competition between the couplings 
in the commensurate and incommensurate antiferromagnetic 
and J- wave superconducting channels. Which of them di- 
verges first may in part depend on the truncation, so when 
one includes fermionic self-energy contributions and quartic 
bosonic couplings, this may have an important effect on the 
size and existence of regions exhibiting local incommensurate 
antiferromagnetic order. 

Most existing studies using the framework of the fermionic 
functional renormalization group focus on the flow of the 
four-fermion vertex, as we do in the present work, so we 
can compare our results with theirs. The most recent re- 
sults for the phase diagram of the two-dimensional Hubbard 
model at varying jj and fixed f, presented in [47], are ob- 
tained by means of an A/^-patch scheme restricted to the flow of 
the four-fermion vertex. A temperature- flow scheme is used 
in that work, where the temperature is used as a parameter 
which flows from infinity to a nonzero value where a first ver- 
tex function reaches some critical value. The temperature T* 
where the divergence of the four-fermion vertex occurs (it is 
obtained in [47] by means of a polynomial fit of the inverse 
susceptibilities) may be compared to our pseudocritical tem- 
perature Tpc. Both correspond, albeit in dififerent ways, to the 
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FIG. 8: Pseudocritical temperatures Tpc for f /t = -0.01 (upper 
panel) and f /t = -0.1 (lower panel) at different values of Ujt - 3.5 
(upper lines), U/t = 3 (middle lines) and U/t = 2.5 (lower lines). 
Short-dashed lines denote the onset of commensurate antiferromag- 
netic order (cAF), long-dashed lines, appearing in small regions at 
larger values of -/u, the onset of incommensurate antiferromagnetic 
order (iAF). Solid lines indicate the onset of J- wave superconducting 
order (dSC)). 

divergence of the four-fermion vertex and the onset of local 
order. However, whereas in the temperature-flow scheme the 
temperature appears as a flow parameter, it is kept fixed in our 
approach. 

The results in pvt] are in complete qualitative agreement 
with ours. This concerns, for instance, the dependence of lo- 
cal order on the values of U and f. If U is increased, the di- 
vergence of vertex functions (equivalent to the emergence of 
local order) occurs already at larger temperatures and the di- 
vergence of the coupling in the J- wave channel, for the small 
values of \t'\ considered, happens closer to the van Hove filling 
[I = Af . There is also agreement on the fact that J- wave su- 
perconductivity, when its coupling is enhanced most strongly, 
competes mostly with incommensurate rather than commen- 
surate antiferromagnetism (cAF) as the dominant instability. 

The quantitative comparison between our results and those 
in ll47ll has to be handled with some care: for T < Tpc our 
flow is stopped at a nonzero scale k. For quantities evaluated 
at ^ 9^ 0, the detailed implementation of the infrared cutofl' 
has an efl'ect on the results. (This contrasts with results for 
temperatures above the pseudocritical line in Fig. [8l where 
the four-fermion vertex never diverges, such that we can ex- 
trapolate to ^ = 0. For ^ = 0, any residual dependence on the 
cutofl' scheme is an indication of the shortcomings of a given 



truncation.) Despite this caveat, the comparison remains in- 
structive. We find that for f It = -0.1 and Ult = 3.5 the 
maximal values of Tpc and T* as functions of jd diff'er by a 
factor of about 4/3, and slightly more for Ujt = 2.5. For 
the onset of J-wave superconducting order, the pseudocrit- 
ical temperature Tpc as a function of ji shown in Fig. [8] is 
larger than T*^^ obtained in ||47|] by a factor of at least 2 for 
both U/t = 2.5 and U/t = 3.5. The diff'erence gets more 
pronounced with increasing distance from half filling. As al- 
ready mentioned, a possible source of quantitative shortcom- 
ings of our calculations is that magnetic fluctuations may be 
overestimated due to our simple parametrization of the mo- 
mentum dependence of the magnetic boson propagator (see 
Eqs. (IB 21) - (IB4I) ). Furthermore, the diagrams shown in Figs. 
[3] and |4] are evaluated at (0,7r) and (tt, 0), which is adequate 
only for not too large values of where the Fermi surface 
is close to the boundary of the Brillouin zone. When mag- 
netic fluctuations are generally overestimated and antiferro- 
magnetic fluctuations are dominant, the critical scales for the 
onset of J-wave superconducting order and hence the pseud- 
ocritical temperatures can be expected to come out too large. 
We recall at this place that both the present work and Ref. 
I47I neglect the renormalization of the fermionic propagator, 
which is expected to have a sizeable lowering efl'ect on the 
value of Tr,c. 



VII. CONCLUSION 

In this work we have shown that the functional renormaliza- 
tion group approach to correlated fermion systems based on 
partial bosonization can account for the competition between 
the antiferromagnetic and superconducting instabilities in the 
two-dimensional Hubbard model. We have studied the emer- 
gence of a coupling in the J- wave channel and its divergence 
in a certain parameter range as a consequence of antiferro- 
magnetic spin fluctuations. In a nutshell, this result confirms 
the spin-fluctuation route to J-wave superconductivity in the 
two-dimensional Hubbard model. 

Our treatment of the fermionic four-point vertex paves the 
way for a unified treatment of spontaneous symmetry break- 
ing in the two-dimensional Hubbard model. In a next step, 
self-energy corrections to the electrons as well as quartic 
bosonic couplings can be included in our approach. This may 
shed light on the question of coexistence of diff'erent types of 
order, which so far has been addressed in the framework of the 
functional renormalization group only on the basis of a mean 
field approach replacing the flow of vertex functions at lower 
scales. 140, llll] In a final step, a unified treatment of the flow of 
vertex functions in both the symmetric and symmetry-broken 
regimes can be given within the present approach. 

Acknowledgments: We are grateful to C. Husemann, A. 
Katanin and M. Salmhofer for useful discussions. SF ac- 
knowledges support by Studienstiftung des Deutschen Volkes. 
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Appendix A: Flowing bosonization 

To illustrate how flowing bosonization works, consider, as 
an example, the part of the efl'ective average action which in 
our truncation involves the m-boson. (The computation is ex- 
actly analogous for the other bosons.) It is given by 

r^ + r^^ + r^= (ai) 

i y m^(-e) {P^(Q) + ml) m(Q) 

Q 



+ Yj h^(K)m(K)-[i/f\Q)ai/f(&)]S(K-Q^&) 

K,Q,Q' 



K\ ,K2,K2i,K4 

Now we introduce a scale dependence of the field m(2), writ- 
ing it as nikiQ). The change in mk(Q) between two scales k 
and k - Ak that are infinitesimally close to each other can be 
chosen to be 



Similarly, we get for the other Yukawa couplings the follow- 
ing contributions from flowing bosonization: 



Pp,k(Q) . 



dkhs(Q) = 5;^/z.[^^.(e)+^r^^^4L,,*(G), (A8) 



^^^dXl AO)- 



Appendix B: Parametrization of bosonic propagators and 
Yukawa couplings 

In our truncation, both the Pi and hi are momentum- 
dependent functions which in addition depend on the scale k. 
In principle, one could try to discretize the momentum depen- 
dence and attempt a numerical solution of the partial difl'eren- 
tial equations for Pi{Q, k) and hi(Q, k). Instead, we proceed in 
this paper to a parametrization of the momentum dependence 
which we describe in this appendix. 



mk(Q) - mk-Ak(Q) = Aak(Q)mk(Q) , 
where the auxiliary field nik(Q) is given by 



(A2) 



(A3) 



and ak(Q) is a function that will be chosen in such a way that 
Ap cancels to zero at all scales. 

To achieve this, we take the generalized flow equation 



(A4) 



which yields in our case 



dk^k = dk^kL + 



Yj {-dkak{Q)Pn,,k(Q)mkiQ) ■ m(Q) 
Q 

+dkak(Q)h„(Q)m,(Q) ■ m(-Q)) .(A5) 



We can read off the modified equations for and and set 
the scale-dependence of to zero, 



dkh„,(Q) = dkK\(Q) + P„,k{Q)dkak(Q), 



(A6) 



This allows us to eliminate the hitherto undetermined function 
Qfy^(2) and to obtain the flow equation for the Yukawa coupling 
including contributions from flowing bosonization. 



dkhm^Q) = ^,/z^L(e)+?P^^,^^l (2). (A7) 



1. Propagators 

The truncation of the inverse bosonic propagators is briefly 
described in the following lines. For a more detailed discus- 
sion see (Note, however, that in [16] we discuss the 
antiferromagnetic propagator which is distinguished from the 
magnetic propagator employed here by a shift in the argument 
by the antiferromagnetic wave vector tt = (tt, tt).) For the ki- 
netic term Pm of the magnetic boson we make the ansatz 



Pm,k(Q) = ^mOjl^A^F(q). 



(Bl) 



The quadratic dependence on frequency is motivated by mean 
field results for small \ojq\. (For larger values of \ojq\, it mim- 
ics the decaying frequency-dependence of the Yukawa cou- 
plings, which is not taken into account explicitly.) 
In Eq. (IBll) we employ for F(q) 



Fc(q) 



(B2) 



if commensurate antiferromagnetic fluctuations dominate. 
Here [q]^ is defined as [q]^ = for qx,y e [-7r,7r] and 

continued periodically otherwise. If incommensurate antifer- 
romagnetic fluctuations dominate, we use 



(B3) 



where the momentum dependence is quartic in momentum 
and explicitly includes the incommensurability q. 



4^2 



F(q, q) = —{(f - [q - irff + 4[q, - nf[qy - nf) . (B4) 
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The shape coefficient used in Eqs. (IB 21) , (IB3I) is com- 
puted as 

= -^(PmiO, 0, 0) - Pm(0, n-q, n)). (B5) 

The Zfn- and -factors are computed from the differences of 
inverse propagators at different frequencies and momenta, 



1 



{Pm,k{27iT, q = 0) - PmM q = 0)) , 



(B6) 



This approximation is most adequate when the loop contribu- 
tions to A'^'^Ki - K2) are minimal for Ki - K2 = and max- 
imal for ^1 = ^2 = n or inversely. This is the case whenever 
either ferromagnetic or commensurate antiferromagnetic fluc- 
tuations dominate. When incommensurate antiferromagnetic 
fluctuations dominate, we approximate the loop contribution 
to hmip(li) occurring in the second line of Eq. (IB 101) by 



where ^ is a parameter which is fixed in such a way that results 
are practically independent of it. For the results displayed be- 
fore, we have set it to 0.15. The propagator of the p-boson is 
treated in exactly the same way. 

For the ^-and J-bosons, the treatment is just as for the a- 
and p-bosons in the commensurate case. Since the minima of 
the inverse propagator do not occur in the vicinity of the wave 
vector TT, we do not use a shift by this vector. 



Ps/d,k(Q) = ^sldCO^ + AsldPsldi^i) 



with 



Psldi^) 



(B7) 



(B8) 



dkhlmm = dkh 



'mlp,k' 



mlpMmtlPk 



(n) 



(BID 



m/p\ 



(n). 



Here Q denotes an incommensurability the size of which is 
given by the positions of the minima of the inverse mag- 
netic (or charge density) propagator. The maximal value of 
hm(0, q), which in the parameter regime we study is always 
located more closely to (tt, tt) than to (0, 0), has been denoted 
/ia in ill 111 H 111]. The minimal value of Pm(Q) corre- 
sponds to the "antiferromagnetic mass term" ml introduced 
there. 



2. Yukawa couplings 

The Yukawa couplings hm(Q) and hp(Q) are parametrized 
by means of a linear momentum dependence 

hm/p(Q) = ^Vr^/^^/pCO) + p^/z^/p(n) , (B9) 

|7r| |7r| 

and the contributions to the flow of hm/p(0) and hm/p(Il) are 
computed according to Eqs. (IA7b and dASb , 

dkhl^u) = ^./^^/p,,L,/,/n)-p./,,,(n)5,4^ 

(BIO) 



For the s- and J-bosons, dominant contributions are for 
Ki = -K3. This is partly accounted for by the propaga- 
tors, and for the ^-boson no further momentum-dependence 
of the Yukawa coupling is assumed. When the J-wave chan- 
nel becomes critical (see Fig. 7 (b) of flT^), the four-fermion 
coupling in the J- wave channel has a sharp peak around zero 
momentum. This is accounted for by including a Gaussian 
function which is centered around zero momentum in the def- 
inition of hd. We have checked that our results are practically 
independent of the width of this Gaussian function, as long as 
it is reasonably peaked. 
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